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We use relativistic Hartree-Fock and correlation potential methods to calculate nuclear spin- 
dependent parity non-conserving amplitudes (dominated by the nuclear anapole moment) between 
hyperfine structure components of the ground state of odd isotopes of K, Rb, Cs, Ba + , Yb + , Tl, Fr, 
and Ra + . The results are to be used for interpretation of current and future measurements. 
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I. INTRODUCTION 

Current best atomic test of the standard model comes 
from the measurements |l| and interpretation |2j of the 
parity-nonconservation (PNC) in atomic cesium (see also 
review 0). The PNC effect measured in cesium is dom- 
inated by the nuclear spin-independent contribution due 
to the weak nuclear charge Qw The value of Q\y( 133 Cs) 
extracted from the measurements is in perfect agreement 
with the standard model. 

Since it is hard to compete with cesium PNC in terms 
of the accuracy of the interpretation of the measure- 
ments, the focus of the atomic PNC studies has shifted 
mostly to the study of the nuclear spin-dependent PNC 
effects (dominated by the nuclear P-odd anapole mo- 
ment) and PNC in the chain of isotopes. The anapole 
moments has been measured so far only for the nucleus 
of 133 Cs [lj]. The information about weak nuclear forces 
extracted from these measurements seems to be inconsis- 
tent with the information obtained from the limit on the 
nuclear anapole extracted from the thallium PNC [U, Q 
and some other sources (see, e.g. |3|). Therefore, it would 
be very important to get more anapole measurements. 
The experiments are in progress or planned for Yb and 
Dy at Berkeley [f| 0] , Yb + at Los Alamos (jl , Ra + at 
KVI HI Q3, Rb and Fr at TRIUMF [TO, [Q. There 
were also proposals to measure atomic PNC in K [3], 
Ba + fl5j . and Xe [l6j]. Most of these experiments are 
aimed at both, nuclear spin-dependent PNC effects and 
PNC in a chain of isotopes. The proposals for K [l4j ]. 
Rb and Fr [ll|, EH specificly target the effect of the 
the nuclear anapole moment in atomic hyperfine transi- 
tions. The interpretation of these measurements requires 
atomic calculations. Note that PNC in the Zeeman tran- 
sitions considered in Ref. [l3[ is also based on such cal- 
culations. 

Nuclear spin-dependent PNC effects in hyperfine tran- 
sitions were considered semiempirically in Refs. (l7l fl8j. 
Ab initio calculations for Fr, which included correlations, 
core pol arization and Breit interaction were reported in 
Ref. [l9j |. Ref. [20[ presents calculations of the nuclear 
spin-dependent PNC for a wide range of single-valence- 
electron atoms. The calculations are in random phase 
approximation (RPA) which means that core polariza- 



tion is included. 

In present paper we perform calculations of the nu- 
clear spin-dependent PNC for a range of atoms of po- 
tential experimental interest. These includes K, Rb, Cs, 
Fr, Tl, Ba+, Ra+, and Yb+. We have included the core 
polarization (RPA) and correlation effects. The results 
in the RPA approximation are in perfect agreement with 
Ref. [2(i| . The inclusion of the correlations in our work 
increases the PNC amplitudes by 3 to 7%. The relatively 
small magnitude of the correlation corrections is due to 
cancelation of different contributions. Based on the de- 
tailed study of the error budget in Ref. 0, we expect 
that the accuracy of present calculations is few per cent. 
The agreement with Ref. for Fr is not so good, and 
we discuss the reasons for the difference. 



II. THEORY 

Hamiltonian describing the nuclear spin dependent 
parity-nonconserving electron-nuclear interaction can be 
written in a form (we use atomic units: h — \e\ = m e — 
1): 



Hi 



SD-PNC 



G F X T , , 



(1) 



where G F w 2.2225 x 10" 14 
of the weak interaction, a = 



is the Dirac ma- 



a.u. is the Fermi constant 
a 
a- 

trie, I is the nuclear spin, and p(r) is the nuclear den- 
sity normalized to 1. The strength of the spin-dependent 
PNC interaction is proportional to the dimensionless con- 
stant k which is to be found from the measurements. 
There are three major contributions to k arising from 
(i) electromagnetic interaction of atomic electrons with 
nuclear anapole moment [21 1 
dependent weak interaction 



(li) electron- nucleus spin- 
|22| , and (iii) combined 
effect of the spin-independent weak interaction and the 
magnetic hyperfine interaction [23| (see also review [|[). 
In this work we do not distinguish between different con- 
tributions to yt and present the results in terms of total 
>c which is the sum of all possible contributions (the nu- 
clear anapole gives the dominating contribution in heavy 
atoms). 
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The PNC amplitude of an electric dipole transition 
between states of the same parity \i) and |/) is equal to: 



El 



PNC 
fi 



E 



(f\d\n)(n\H PNC \i} 



Ei — E n 
(f\H PNC \n)(n\d q \z) 



E f ~E n 



(2) 



where d — — e J2i r i is the electric dipole operator, \a) = 
| J a F a M a ) and F = I+J is the total angular momentum. 

Applying the Wigner-Eckart theorem we can express 
the amplitudes via reduced matrix elements 



E1% NC 



= (-l) F f- M f 
x (JfFf 



F f 1 Fi 
-M f q Mi 

dpNc\\JiFi)- 



(3) 



Detailed expressions for the reduced matrix elements of 
the SI and SD PNC amplitudes can be found e.g. in 
Refs. [HI and [ij. For the SD PNC amplitude we have 

(Jf,F f \\d BO \\J i ,Fi) = ^X 

xyj(l+ 1)(2J + l)(2Fi + l)(2F f + I)/ 7 



e [(-d*-* { t J k 

(J f \\d\\n, J n )(n, J n \\ap\\Ji) 

E n — Ei 

(J f \\ap\\n, J n )(n, J n \\d\\Ji) 



J n Jf 1 
F f Fi I 



(4) 





Jf 




/ Jn 
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Ff } 




Ff I ] 



E n - Ef 



This formula can be used for optical and microwave tran- 
sitions. For the microwave transitions it can be further 
simplified as it has been done in Ref. [l9| . 

To calculate the PNC amplitudes we use direct sum- 
mation over a complete set of single-electron states con- 
structed using the B-spline technique [24j. Correlations 
and core polarization effects are included by means of the 
correlation potential method (25|. The calculations are 
very similar to those of our previous work 26] , therefore 
we omit here the details. 



III. RESULTS 

The results of calculations are presented in Table Q] 
The results of present work for the PNC amplitude are 
presented in two different approximations, the RPA ap- 
proximation, which includes core polarization but no cor- 
relations beyond it, and final results which include core 
polarization and Brueckner-type correlations. The re- 
sults of Ref. [2(|, which were obtained in the RPA ap- 
proximation, are also presented for comparison. For the 



convenience of the comparison we present the results for 
the PNC amplitudes in a form which corresponds to the 
PNC Hamiltonian, used in (20| : 



#sd-pnc = -j=talp{r) 



(5) 



It differs from (p} by the use of the different weak inter- 
action constant £ (£ = a 1 1) ■ 

The RPA results of 20] and present work are in perfect 
agreement with each other. The inclusion of Brueckner- 
type correlations increase the PNC amplitudes by 3 to 
7%. It is interesting to note that the correlation correc- 
tion is larger for light elements and slowly decrease for 
higher Z. 

In Table Q] we also present the values of the reduced 
matrix elements for the magnetic dipole (Ml) transition 
amplitudes between corresponding hyperfine states and 
the degrees of circular polarization of light P. The Ml 
amplitudes are given by 



(F,J\\M1\\F',J) = (-l) I+J + F g(J,L)p B (6) 

F 1 F' 
J 



: V(2F + 1)(2F' + 1) J( J + 1)(2J + 1) | F j J 



where g( J, L) is the ^-factor of atomic state J,L,S= 1/2 
J(J+l)-L(L + l) + § 



g{J,L) = l + 



2J(J+ 1) 



(7) 



J, L are atomic total and angular momentums, and \ib 
\e\H/2mc is Bohr magneton. In the case of J = 1/2, F 
I ±1/2 the Eq. H can be further reduced to 



Ml 



1 1(1 + 1) 



L(L+1) 



(8) 



The degree of the circular polarization of light is given 
by 



P = 2 



Im(El) 
Ml 



(9) 



The results for 211 Fr, converted to a different defini- 
tion of the weak interaction constant {x = £/) consis- 
tent with Hamiltonian (JTJ, are 0.529 x lQ~ w ix H?J and 
0.553 x 10^ 10 ix (present work). They differ by only 
4.5% and this difference is due to correlations which 
were included in the present work but not included in 
Ref. [13] ■ On the other hand, the difference between the 
result of present work and the calculations by Porsev 
and Kozlov [19(, which is 0.491 x lO^ 10 ^, is significantly 
larger, being about 13%. It is important to understand 
the reason for this difference, since the calculations by 
Porsev and Kozlov are the only other calculations for 
Fr which included correlations beyond the RPA approx- 
imation. The experimental work for Fr is in progress at 
TRIUMF [TO, Oj] and for its future interpretation it is 
important to have reliable theoretical results. 
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TABLE I: Reduced matrix elements of the nuclear spin-dependent weak interaction for microwave transitions between ground- 
state hyperfine levels F l ,F2 in atoms and ions of potential experimental interest. In this table, Z is the nuclear charge, A is 
the atomic number, I is the nuclear spin, the hyperfine levels have angular momentum F = I ± 1/2, Mi is magnetic dipole 
transition amplitude, P is the degree of circular polarization. All dimensional numbers are in atomic units. The values for 
Ml amplitudes include Bohr magneton (/j,b ~ a/2 in atomic units). The PNC amplitudes and circular polarization P are 
proportional to the weak interaction constant £ (0) which is omitted from the table. Numbers in square brackets represent 
powers of 10. 



Atom 


Z 


A 


I 


Fi 


F 2 


RPA 
Ref. [20] 


RPA 


RPA+E 

this work 


Mi 


P 


K 


19 


39 


3/2 


2 


1 


2.222[-13] 


2.222[-13] 


2.378[-13] 


0.0100 


4.76[-ll] 


Rb 


37 


85 


5/2 


3 


2 


2.550[-12] 


2.532[-12] 


2.660[-12] 


0.0125 


4.27[-10] 


Rb 


37 


87 


3/2 


2 


1 


1.363 [-12] 


1.353[-12] 


1.422[-12] 


0.0100 


2.85 [-10] 


Cs 


55 


133 


7/2 


4 


3 


1.724[-11] 


1.720[-11] 


1.819[-11] 


0.0145 


2.51 [-09] 


Ba+ 


56 


135 


3/2 


2 


1 


6.169[-12] 


6.186[-12] 


6.385[-12] 


0.0100 


1.28 [-09] 


Yb+ 


70 


171 


1/2 


1 







3.844[-12] 


3.987[-12] 


0.0063 


1.26 [-09] 


Yb+ 


70 


173 


5/2 


3 


2 




2.274[-ll] 


2.359 [-11] 


0.0125 


3. 79 [-09] 


Tl 


81 


203 


1/2 


1 





3.000[-ll] 


3.001[-11] 


3.066[-ll] 


0.0021 


2. 92 [-08] 


Fr 


87 


211 


9/2 


5 


4 


2.379[-10] 


2.362[-10] 


2.487[-10] 


0.0162 


3.07[-08] 


Fr 


87 


223 


3/2 


2 


1 


5.820[-ll] 


5.815[-11] 


6.051[-11] 


0.0100 


1.21 [-08] 


Ra+ 


88 


223 


3/2 


2 


1 


5.987[-ll] 


5.991[-11] 


6.147[-11] 


0.0100 


1.23 [-08] 



TABLE II: Contributions to the 6s — 7s -Epnc for Cs in units 
10 _11 iea.B(— Qw/N). First line is equivalent to the approxi- 
mation used in present work. 



Brueckner orbitals + core polarization 0.9077 

Weak correlation potential 0.0038 

Structural radiation 0.0029 

Normalization -0.0066 

Breit correction -0.0055 

Neutron distribution -0.0018 

Radiative corrections" -0.0029 

Total 0.8976 



Ref. [33 



The most obvious difference between present calcula- 
tions and those of Ref. is inclusion of higher-order 
correlations in our work and some small effects in [l9j . 
These small effects include Breit interaction, structural 
radiation and renormalization of wave functions. We first 
discuss these small corrections. The most detailed study 
of all important corrections to the PNC amplitude has 
been done for the 6s - 7s PNC amplitude for Cs @,IH- 
The results are summarized in Table [TXJ The relative val- 
ues of specific corrections for the hfs PNC transition in 
Fr might be slightly different, however, the qualitative 
picture should be very similar. First line in Table |TT] cor- 
responds to the approximation used in the present work. 
Note that the final result is only about 1% smaller. Fur- 
thermore, there is strong cancelation between different 
contributions. For example, the contributions from the 
weak correlation potential and the structural radiation 
are canceled by the renormalization of the wave func- 
tions. Different contributions to the radiative corrections 
are not so small (up to -0.8%). However, they have dif- 
ferent signs and nearly cancel each other. This means 



TABLE III: Reduced matrix elements of the nuclear spin- 
dependent PNC amplitude (7s, F = 5|jd PN c||7s, F = 4) for 
211 Fr in units xl0" 10 ixr. To compare these numbers with 
those presented in Table U one should multiply them by the 
nuclear spin / = 9/2. Contributions from the core are not 
included in first three columns. 





DHF 


RPA 


MBPT" 


core 


Total 


Porsev and Kozlov [19] 
this work 


0.418 
0.423 


0.058 
0.058 


0.033 
0.009 


-0.018 
0.063 


0.491 
0.553 



"Many-body perturbation theory for the correlation correction. 



that inclusion of some of the small corrections while not 
including others cannot be justified and may lead to less 
accurate results. 

We proceed to the examining the correlations. For this 
purpose we compare our results with those from Ref. [ilj ] 
term by term as it is shown in Table [TTTT The first column 
shows the PNC amplitude in the Dirac-Hartree-Fock ap- 
proximation without any correlations. There is a small, 
about 1% difference in the results. Some of this differ- 
ence might be attributed to the Breit interaction which 
is included in [l9| and not included in our work. There 
might be also some difference due to different treatment 
of the nucleus. We use smooth Fermi distribution of the 
nuclear electric charge and anapole moment, while a step- 
like function is used in (l9j . 

The RPA correction is the same in both works (see 
Table IIIII) . The difference in correlation correction is 
significant which is not surprising due to very different 
treatment of the correlations in two works. We include 
the dominating Brueckner-type correlations with the use 
of the all-order correlation potential E (2(| . The correla- 
tions are included in Ref. [19( in the second-order only, 
including the structure radiation and the renormalization 
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of the wave functions. These latter corrections are small 
(see discussion above) and cannot explain the difference 
in the results. 

To test whether the difference is due to the higher- 
order correlations we performed calculations in which 
all higher-order correlations were removed. The corre- 
lation potential £ was calculated in the second order 
and the core polarization corrections were not included. 
The result, 0.021 x 10~ 10 i>c, is in a reasonable agreement 
with [l|| . In the end, the difference in the results due to 
different treatment of the correlations is about 5%. This 
difference is most likely due to the higher-order correla- 
tions included in the present work. 

The largest difference comes from the contribution of 
the core states. This contribution is small and negative 
in Ref. [l9|. In our calculations it is not so small and it 
is pos itive (see Table llllj) . The value of this correction 
in [l9| suggests that RPA corrections were probably not 
included. Our value for the core contributions without 
RPA corrections is +0.015 x 10~ 10 ix. However, there is 
no clear explanation for the different sign. Note, that we 
do not distinguish between core and excited states in the 
summation over complete set of states in This leaves 
no room for a sign error. On the other hand, to the best 
of our knowledge, the contribution of the core states was 
calculated separately in Ref. [r| . 

There is a simple test to check the sign of the core 
contribution. The summation over core states in (j4]) is 
dominated by the &P1/2 state, while summation above 
core is dominated by the 7pi/ 2 state. The angular coeffi- 
cients in ((U) are the same for core and higher states. One 
needs only to compare the energy denominators and the 
matrix elements of weak and electric dipole interactions. 
The energy denominators E^ s — Eq Pi/2 and Ej s — Ef Pl/2 
are of the opposite sign. To compare the signs of matrix 
elements we need to fix the phase of the wave functions. 
It is convenient to have f(r) > at r — > 0, where f(r) 
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is the upper component of the Dirac spinor. Then, the 
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mation. The many body corrections are not large enough 
to change the sign of the electric dipole matrix elements. 

In our previous work [28|, we used the results of 
Ref. [IH and the ratio of the matrix elements of spin- 
dependent PNC interaction and electron electric dipole 
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hancement factor for Fr from the spin-dependent PNC 
calculations of 19]. The result, d(Fr) = 854d e , was in 
reasonably good agreement with the many-body calcu- 
lations of Ref. 29], d(Fr) = 910(46]d e . We can use the 
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